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For any translation-invariant quantum lattice system with a symmetry group G, we propose a practical and 
universal construction of order parameters which identify quantum phase transitions with symmetry-breaking 
order. They are defined in terms of the fidelity between a ground state and its symmetry-transformed counterpart, 
and are computed through tensor network representations of the ground-state wavefunction. To illustrate our 
scheme, we consider three quantum systems on an infinite lattice in one spatial dimension, namely, the quantum 
Ising model in a transverse magnetic field, the quantum spin- 1/2 XYX model in an external magnetic field, 
and the quantum spin-1 XXZ model with single-ion anisotropy. All these models have symmetry group Z 2 
and exhibit broken-symmetry phases. We also discuss the role of the order parameters in identifying factorized 
states. 
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Introduction. In the conventional Landau-Ginzburg-Wilson 
paradigm a basic notion is spontaneous symmetry breaking, 
which is traditionally characterized in terms of a local order 
parameter |l|, |2fl. Typically, such a local order parameter is 
model-dependent and not always obvious to define. Although 
attempts have been made to derive a local order parameter 
for quantum lattice systems undergoing quantum phase transi- 
tions (QPTs) JslH], it is highly desirable to find a simple, uni- 
fying, and model-independent way to characterize symmetry- 
breaking order when it occurs in a quantum lattice many -body 
system. 

In this Letter, we address this issue from a quantum in- 
formation perspective through the notion of fidelity. In 
Refs. J5I H HUh it has been argued that the ground-state 
fidelity per site may be used to detect QPTs. Since this ar- 
gument is based solely on a basic postulate of quantum me- 
chanics regarding quantum measurement, the approach is ap- 
plicable to quantum lattice systems in any number of spatial 
dimensions, regardless of the type of internal order present in 
quantum many-body states. It has been confirmed that the 
ground-state fidelity per site is able to identify QPTs aris- 
ing from spontaneous symmetry breaking (J, la LZL IM |9|]> the 
Kosterlitz-Thouless transition [10], and topological QPTs in 
the Kitaev model ll ill . In an extension of this notion, here we 
propose a universal approach to define and compute order pa- 
rameters for any translation-invariant quantum lattice system, 
with a symmetry group G, undergoing a QPT with symmetry- 
breaking order. We perform the explicit computation of order 
parameters through the use of tensor network al gori thms, in 
particular the matrix product states (MPS) HHHS [3 for 
systems in one spatial dimension. 

To illustrate our scheme, we investigate the following mod- 
els on an infinite lattice in one spatial dimension: the quantum 
Ising model in a transverse magnetic field, the quantum spin- 
1/2 XYX model in an external magnetic field, and the quan- 
tum spin-1 XXZ model with single-ion anisotropy. All these 
systems possess a discrete symmetry group Z2, which is spon- 



taneously broken as the system undergoes a QPT. Although 
these examples are restricted to Z2-systems on an infinite-size 
lattice in one spatial dimension, we emphasize that the scheme 
extends to any translation-invariant quantum lattice system, 
with arbitrary symmetry group, in any spatial dimension. For 
systems in higher spatial dimensions the computation of the 
ground state is accommodated by the tensor product states 
(TPS) [ 15], or equivalently, the projected entangled-pair states 

(PEPSMlffl. 

The construction of these order parameters not only allows 
us to locate critical points, but also enables us to identify fac- 
torized states \t//(Af)), where Af is the so-called factorizing 
field 118, [l9ll . Such a factorized state FF(A/)) can occur in 
the symmetry-broken phase. The fact that no entanglement 
exists makes it the most ordered state, with a salient feature 
that the order parameters take their maximum value at A — Af. 

Universal construction of order parameters for translation- 
invariant quantum lattice systems with symmetry-breaking or- 
der. Consider an infinite-lattice translation-invariant quantum 
system with symmetry group G and Hamiltonian H(A), with A 
a control parameter. According to Wigner's theorem the rep- 
resentations of g e G are either unitary or anti-unitary. The 
group elements which are represented unitarily form a sub- 
group G c G, and we hereafter restrict our considerations to 
this subgroup. Suppose a symmetry-breaking QPT occurs at a 
critical point A c where, without loss of generality, we assume 
there is no symmetry breaking for A > A c . For the ground 
state \tf/(A)) in the symmetric phase, the fidelity \(if/{A)\g\ip(A))\ 
is equal to one for any symmetry operation g e G. In the 
broken-symmetry phase A < A c we can have 



< wx>WW)\ 2 



(i) 



where \iJ/(A)) and \if/'(A)} are any two states in the degenerate 
ground state. This description is valid for any system admit- 
ting a QPT with symmetry-breaking order, regardless of the 
type of symmtery group or whether the transitions are con- 
tinuous or discontinuous. For later use we will denote the 
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FIG. 1 : (color online) The computation of the transfer matrix E from 
the infinite Matrix Product State (iMPS) representation of ground 
state wavefunctions for translation-invariant quantum system on an 
infinite one-dimensional lattice. The transfer matrix E is constructed 
from three-index tensors Ta, T b , T' a , and T' B attached to sites, and di- 
agonal singular-value matrices A A , A B , A' A , and A' B attached to bonds, 
together with two square gates representing a nontrivial element g of 
the symmetry group G acting on physical indices. 

ground-state subspace as V(A). 

Based on this observation we now describe a procedure to 
construct order parameters which provide quantitative infor- 
mation about symmetry-breaking QPTs. Because QPTs gen- 
erally occur in the infinite lattice limit, we will utilize an al- 
gorithm developed by Vidal [17]. This algorithm provides 
an efficient way to generate an infinite MPS (iMPS) repre- 
sentation of the system's ground-state wavefunction, and in 
turn compute expectation values. Determination of the ground 
state amounts to computing the imaginary time evolution op- 
erator exp(-i/r) acting on an initial state 1^(0)): ^(t)) = 
exp(-i/r)| v P(0))/|exp(-//T)| x P(0))|. In the symmetry-broken 
phase with ground-state degeneracy, we can obtain more than 
one ground-state representative through different choices of 
the initial state PPfTj)). Let \i//(A)) and \ip'{A)) denote two such 
iMPS representations (which may or may not be equivalent 
states). The structure associated to the algorithm, which ut- 
lizes the translational invariance of the system, leads to the 
conclusion that 

mX)\g\il/(X))\ = lim |tr(£ L )|, (2) 

L— >oo 

where E, referred to as the transfer matrix, is the four-index 
tensor schematically defined in Fig. [TJ in terms of diago- 
nal singular-value matrices Aa, Ab and the three-index tensors 
Ia, Tg. The iMPS representation becomes exact in the limit 
as the dimension of the matrices Aa, Ab approaches infinity. In 
practice, we introduce a truncation dimension^ associated to 
Aa,Ab and adjust^ to achieve an extrapolation of results for 
the^ — > oo case. 

Observe next that (f2]) can only take the values zero or one, 
since the trace is simply the sum of the eigenvalues. This 
is in stark contrast to (TTJ, and indicates that the iMPS algo- 
rithm is preferential in its convergence to states in the ground- 
state subspace. Fixing a representative \if/j(A)), we define 
Wj(A) = span{\i[/f,j(A)) = h\t//j(A)) : h e G}. Remembering 



that the representation of G is unitary, it also follows that the 
representation of each g e G on Wj(A) is of the form of a 
permutation matrix (up to a diagonal unitary transformation). 
The existence of such a representation for all groups G fol- 
lows from Cayley's theorem, which states that every group 
is the subgroup of a symmetric group. It may be that Wj{A) 
does not span the degenerate subspace of the broken symme- 
try phase. Within the iMPS algorithm, we can generate sub- 
spaces orthogonal to Wj(A) simply by choosing the initial state 
[¥(0)) of the algorithm to be orthogonal to Wj(A). In principle 
a basis for V(A) can be constructed in this manner leading to 
V(A) = ®jWj(A). We then still have that the representation of 
geGon V(A), which has a block diagonal structure, is of the 
form of a permutation matrix. 

We can now exploit these facts to uniquely define order pa- 
rameters which characterize the broken-symmetry phase. For 
a fixed choice of g e G and a fixed iMPS representation 
\if/(A)), let f g (A) denote the square root of the eigenvalue of 
E which has the largest absolute value. As such, one sees that 
|/#(/i)| = 1 for all g e G in the symmetric phase A > A c , and 
< |/ g (/l)| < 1 in the broken-symmetry phase A < A c . We 
now define I g (A) to be 

I g (A) = ^1 - \f g (A)\ 2 . (3) 

and argue that the set O = {I g (A) : g e G} defines a set of 
order parameters which detect QPTs with symmetry -breaking 
order. First, each I g (A) is zero if A > A c . Second, each I g {A) 
can take a value ranging from to 1 if A < A c . These fea- 
tures are nothing but what we require for I g (A) to be an order 
parameter. Moreover, O = {I g (A) : g e G) is unique on each 
block representation Wj(A) c V{A), in that the set is indepen- 
dent of the choice of iMPS representation \\j/(X)) e Wj(A). To 
show this, we first consider for each G-congugacy class C the 
set of order parameters Oc - {I g (A) : g e C). Any different 
iMPS representation \ip'j(A)) e Wj(A) will be related to \if/j(A)) 
by \4f'j(A)) - exp(i(p)\if/h,j(A)) for some h e G and some phase 
<p, due to orthogonality of \iffj(A)) and \ift'j(A)) as decreed by 
(O. The uniqueness of the set Oc for Wj(A) follows from 

\(i//j(A)Wj(A))\ = mj(A-)Wh,M))\ = \(<fij(A)\h-'gh\<J,j(A))\. 

Consequently, O = UcOc is a uniquely-defined set of order 
parameters on each subspace block Wj(A) c V(A). This is in 
some contrast to the conventional notion of an order parame- 
ter, whose value is associated with individual states in V(A). 
These considerations are valid for any quantum lattice system 
where there is symmetry-breaking order, thus the construction 
can be applied universally. Finally, we define a symmetry- 
breaking QPT to be second-order if all I g (A) are continuous 
functions of A. If there are any discontinuous I g (A) we say the 
QPT is first-order. 

The models. To illustrate our scheme, we consider three 
one-dimensional lattice Hamiltonians with Z2-symmetry. 
This group is generated by a single non-trivial element g, 
which squares to the identity, and in each case admits a unitary 
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FIG. 2: (color online) The universal order parameter 1(A) for the 
quantum Ising model in a transverse magnetic field. For the con- 
trol parameter A less than the critical value A c = I, the Z2 symmetry 
is spontaneously broken. In the iMPS simulation, a transition point 
A x occurs for a given truncation dimension x- As the truncation di- 
mension x is increased, A x approaches the exact value. Note that 
the factorized field Af occurs at Af = 0, at which the universal or- 
der parameter takes the maximum value. Inset: the critical point A c 
is determined from an extrapolation of the pseudo phase transition 
point A x with respect to the truncation dimension x- Here, the fit- 
ting function is A x = Ac + ax~ b , with A c = 1.0023, a = 0.3937 and 
b = 2.4390. 
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FIG. 3: (color online) The universal order parameter 1(A) for the 
quantum XY model in a transverse magnetic field. For increasing 
A, a transition point A x occurs as 1(A) goes to zero. With increasing 
truncation dimension x, the transition point A x approaches the crit- 
ical point Ac = 1. The universal order parameter 1(A) reaches the 
maximum value when Af is equal to 0.5. Inset: the critical point A c 
is determined from an extrapolation of the pseudo phase transition 
point A x with respect to the truncation dimension x- Here, the fitting 
function is A x = A c + ax~ b , with A c = 1.000039, a = 0.197301 and 
b = 1.845012. 



representation. The construction described above thus leads to 
a single order parameter, which we hereafter denote as 1(A). 

The first of the models is the quantum XY model in a trans- 
verse magnetic field in an infinite-size lattice in one spatial 
dimension. The Hamiltonian takes the form, 



1 + y 



S®S 



[i+i] 



1-7. 



(4) 



where S„ (a = x,y,z) are the spin- 1/2 Pauli operators at 
lattice site i, A is a transverse magnetic field, and y is an 
anisotropic coupling constant. The model is invariant under 
the symmetry operation: S% — > —S^\ Sy — » —S y and 
— > 5j ] for all sites, which yields the Z2 symmetry. For 
nonzero y, it is critical at A = 1 1I20I1 . In addition, a factor- 
izing field occurs at Af = yl - y 2 . If y - 1, the quantum 
XY model reduces to the quantum Ising model in a transverse 
field. 

We also investigate the quantum spin- 1/2 XYX model in an 
external magnetic field. The Hamiltonian can be written as 



H = J] (S®SW + A yS fs\! +1] + SfS^ 1] + hS^ ] ) , (5) 
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FIG. 4: (color online) The universal order parameter 1(h) for the spin- 
1/2 quantum XYX model in an external magnetic field. A transition 
point h x occurs as as 1(h) vanishes. As the truncation dimension 
X increases, the transition point h x approaches the critical point h c . 
In addition, the universal order parameter 1(h) takes the maximum 
value at hf ~ 3.16. Inset: the critical point h c is determined from 
an extrapolation of the pseudo phase transition points h x with respect 
to the truncation dimension x- Here, the fitting function is h x = 
h c + ax~ b , with h c = 3.2052, a = 0.0277 and b = 0.9848. 



where S„ (a = x,y,z) are the Pauli spin operators at site 
i, A y is a parameter describing the rotational anisotropy, and 
h is an external magnetic field. This model also possesses 
a Z2 symmetry, with the symmetry operation: — > — , 



and SJ? 



'I 



S 1 , for all sites. Below we shall 



choose A v = 0.25. In this case, the critical magnetic field is 



h c ~ 3.210(6) [ 18], with a factorizing field h f ~ 3.162. 



The last model considered here is the quantum spin 1 XXZ 
model with single-ion anisotropy. The Hamiltonian takes the 
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FIG. 5: (color online) The universal order parameter 1(D) for the 
quantum spin-1 XXZ model with single-ion anisotropy, which ex- 
hibits discontinuity at the transition point D x . This indicates that 
a first-order QPT occurs at D = D x . As the truncation dimen- 
sion^ increases, the transition point D x approaches the critical point 
D c = 9.9434. Inset: the fitting function is D x = D c + ax' b , with 
n = 9.9434, a = 1.1490 and b = 1.4033. 



form, 

oo oo 

H = J] [J(S 1 M + 1] + SfS^) + J z SfS^ +1] ]+D J] Sf 2 , 

i— — CO /= — CO 

(6) 

where (« = x,y,z) are the spin-1 operators at the i-th 
lattice site , D represents uniaxial single-ion anisotropy. We 
choose / = 1, J z = 10, with D as the control parameter. As 
such, the system undergoes a first-order QPT from a gapped 
Z2 symmetry-broken Neel phase to a gapped Z2 symmetric 



i] 



large-D phase, with the symmetry operation: S 
Sf -Sf and S® -S® for all sites. 

The results. In Figf2] we plot the universal order parameter 
1(A) for the quantum Ising model in a transverse field, with 
the field strength A as the control parameter. For the control 
parameter A less than a pseudo critical value A x , the univer- 
sal order parameter 1(A) is non-zero, which characterizes the 
Z2 symmetry-broken phase. In the symmetric phase A > A x , 
the universal order parameter 1(A) is zero. When the control 
parameter A varies across the point A x , the behavior of the uni- 
versal order parameter 1(A) implies that the system undergoes 
a phase transition at the A x . As the truncation dimension^ is 
increased, A x moves toward the known critical point A c = 1 . 
Performing an extrapolation of A x with respect to^-, we obtain 
A c = 1.0023. Note that the universal order parameter reaches 
the maximum value at A = 0. The maximum value of 1(A) 
coincides with the existence of a factorized state at this point, 
in which no entanglement exists. 

In Figf3] the universal order parameter 1(A) for the quan- 
tum spin 1/2 XY model in an external magnetic field A is plot- 
ted, with the magnetic field strength A as the control param- 
eter. A transition point A x occurs as 1(A) varies from zero to 
nonzero values. With increasing trucation dimension the tran- 
sition point A x approaches the critical point A c - I. Perform- 



ing an extrapolation with respect to x yields the critical point 
A c - 1.000039. A factorizing field A f occurs at A f = 0.5, 
where again 1(A) takes its maximum value. 

In Fig|U we show the universal order parameter 1(h) be- 
tween a ground state and its symmetry-transformed counter- 
part for the quantum XYX model in an external magnetic 
field. In the range h < h x , the universal order parameter 1(h) is 
non-zero, which characterizes the Z2 symmetry-broken phase. 
A transition point h x occurs as 1(h) vanishes. As the trunca- 
tion dimension ;f increases, the transition point h x approaches 
the critical point h c = 3.2052. Performing an extrapolation 
of h x with respect to^, we obtain h c - 3.2052, a = 0.0277, 
b = 0.9848. The universal order parameter 1(h) takes the max- 
imum value at hf ~ 3.16. 

In Figj5] the universal order parameter 1(D) is plotted for 
the quantum XXZ model with single-ion anisotropy. For D < 
D x , the system is in the gapped Z2 symmetry-broken Neel 
phase, so the universal order parameter 1(D) is nonzero. For 
D > D x , it is in the gapped Z2 symmetric large-D phase, so 
1(D) is zero. We see that 1(D) is discontinuous, providing an 
example of a first-order QPT. As the truncation dimension x 
increases, the transition point D x approaches the critical point 
D c = 9.9434. 

Summary. We have described a universal procedure to com- 
pute order parameters for any translation-invariant quantum 
lattice system with a symmetry group G. The scheme has 
been illustrated for the quantum Ising model in a transverse 
magnetic field, the quantum spin- 1/2 XYX model in an ex- 
ternal magnetic field, and the quantum spin-1 XXZ model 
with single-ion anisotropy, by exploiting the infinite MPS al- 
gorithm for one-dimensional systems. In all instances the pro- 
cedure is successful in identifying the QPT points. Moreover, 
we observe that occurences of factorizing fields are in corre- 
spondence with the maximal values of the order parameters. 
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